Properties of the superconducting state in compressed Sulphur 
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The thermodynamic properties of the superconducting state in Sulphur under the pressure at 160 
GPa were determined. It has been shown that: (i) the critical value of the Coulomb pseudopotential 
is equal to 0.127; (ii) the critical temperature {Tc = 17 K) should be calculated by using the modified 
Allen-Dynes formula; (iii) the effective electron-electron interaction is attractive in the range of 
frequencies from zero to the frequency slightly lesser than the maximum phonon frequency (~ 
CSSfimax); (iv) the dimensionless ratios 2A (0) /fesTc, AC (Tc) /C^ (Tc) and TcC^ (Tc) /Hi (0) 
are equal to 3.7, 1.65 and 0.16 respectively; (v) the ratio of the effective to bare electron mass 
reaches maximum of 1.77 for T = Tc- 

PACS numbers: 74.20. Fg, 74.25. Bt, 74.62. Fj 



I. INTRODUCTION 

At ambient pressure (p) Sulphur (S*) is the insulator. 
Under compression S undergoes the sequence of struc- 
tural phase transitions and metallizes [1|. At 93 GPa, 
Sulphur is the superconductor with Tc = 10.1 K 2]. 
The critical temperature increases linearly with the pres- 
sure up to 157 GPa; the rate of dTc/dp is equal to 0.055 
K/GPa and [Tc] igygp^ = 14 K. Near 160 GPa, when 
5" transforms to the /3-Po-type phase, the critical temper- 
ature increases rapidly from 14 K to 17 K. We note that 
Tc = 17 K is among the highest critical temperatures ob- 
served in the elemental solids (the maximum value of Tc, 
which is equal to 25 K, has Calcium under the pressure 
at 161 GPa [2|]). 

In the paper we have studied the thermodynamic prop- 
erties of the superconducting state in S for p — 160 GPa. 
The numerical calculations have been conducted in the 
framework of the Eliashberg formalism Q . 



II. THE ELIASHBERG EQUATIONS 

The Eliashberg set in the mixed representation takes 
the form: 
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FIG. 1: The critical temperature as a function of /i*. The 
results have been obtained with an use of the following ap- 
proaches: the Eliashberg equations, the modified (classical) 
Allen-Dynes formula and the McMillan expression. The ar- 
row shows the experimental value of Tc for ^c = 0.127. 
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FIG. 2: (A)-(B) The order parameter and the wave function 
renormahzation factor on the imaginary axis for the selected 
temperatures. (C)-(D) The dependence of Am=i and Zm^i 
on the temperature. The values of Am=i (T) can be fitted by 



the expression: Am=i (T) — Am= 
A„=i (0) = 2.7 meV and /3 = 3.37. 
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where 



The symbols (j) i'^) and (pm = (j) {i^m) denote the or- 
der parameter functions on the real and imaginary axis 
respectively; Z (w) and Zm = Z {iuJm) are the wave func- 
tion renormalization factors; m-th Matsubara frequenc;^ 
is given by: uj„i = {n/(3) {2m — 1), where /3 = [ksTY 
and kB is the Boltzmann constant. In the framework 
of the Eliashberg formalism, the order parameter is de- 
fined as: A = (p/Z . The electron-phonon pairing kernel 



has the form: \{z) = 2 J^^ 



dQ-, 



^-F (il), where 



the Eliashberg function for Sulphur under the pressure 
at 160 GPa {a^F [Vl)) has been calculated in the paper 
0; the maximum phonon frequency (fimax) is equal to 
86.7 meV. The function /j,* (a;™) = ^*9{ujc — \uJm\) de- 
scribes the electron depairing interaction; 9 denotes the 
Heaviside unit function and Wc is the cut-off frequency 
(ujc = 3riniax)- The critical value of the Coulomb pseu- 
dopotential (/i^) should be calculated by using the con- 
dition: Am=i {^J.*) =0 for T = Tc. The following result 
has been obtained: /i^ = 0.127. The symbols N {cu) 
and / (cj) denote the statistical functions of bosons and 
fcrmions respectively. 

The Eliashberg equations have been solved for 2201 
Matsubara frequencies (M = 1100) by using the numeri- 
cal method presented in the paper [6] . In the considered 
case the functions 6 and Z are stable for T > 4 K. 



III. RESULTS 

The classical Allen-Dynes or McMillan formula cannot 
be used in the case of Sulphur (see Fig. [I]) [!]■ For 
this reason, we have obtained modified Allen-Dynes ex- 



FIG. 3: The dependence of the real and imaginary part of the 
order parameter on the frequency for selected temperatures. 
The rescaled Eliashberg function is also plotted. 




FIG. 4: The order parameter on the complex plane for the se- 
lected values of the temperature. The bold lines represent the 
solutions for uj £ (0, fimax); whereas regular lines correspond 
to the solutions for u) £ (f2max,<^c). 



pression on the basis of 330 exact values of Tc (m*)- Iii 
particular: 
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where fi is the strong-coupling correction function (/i 
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) and /2 denotes the shape correction func- 



The symbol uj2 denotes the 



second moment of the normalized weight function, uiin 
is the logarithmic phonon frequency and A is called the 
electron-phonon coupling constant. For Sulphur the pa- 
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FIG. 5; (A) The dependence of the specific heat in the su- 
perconducting and normal state on the temperature. (B) The 
thermodynamic critical field as a function of the temperature. 



rameters -^o^, i^in and A are equal to 45.39 meV, 37.7 
meV and 0.75 respectively. The functions Ai and A2 

have the form: Ai = 55/x* - 1 and A2 = 11.06/Li* . 

In Figs. [5] (A) and (B) we have shown the solutions 
of the Eliashberg equations on the imaginary axis for 
the temperature range from 4 K to 17 K. Additionally, 
in Figs. [2] (C) and (D) the functions A„=i (T) and 
Zm=i (T) are plotted. On the basis of the presented re- 
sults one can state that the value of the order param- 
eter at the temperature of zero Kelvin (A(0)) is equal 
to A (T = 4K) with the good approximation. In order 
to calculate A (T = 4K) the following algebraic equation 
has to be used: A (T) = Re [A (w = A (T))]. The form 
of the order parameter on the real axis for T = 4 K is 
presented in Fig. [3] (A). With the help of the simple cal- 
culation one can obtain: A{T — 4K) = 2.71 meV. Then, 
the value of the ratio Ri = is equal to 3.7. We 

notice, that in the framework oftfie BCS model, the pa- 
rameter Ri takes the lower value: [-RiIbcs ~ 

3.53 i. 

On the basis of Figs. [3] (A)-(D) we have stated that 
Re[A (w)] and Im[A (w)] are plainly correlated with the 
shape of the Eliashberg function in the full range of the 
considered temperatures. From the physical point of view 



this indicates, that the form of order parameter on the 
real axis clearly reflects the form of the electron-phonon 
interaction in Sulphur. Additionally, in Fig. |4] we have 
plotted the order parameter on the complex plane for the 
selected temperatures; the frequencies from to uJc have 
been taken into consideration. We have found that the 
values of A (w) form the distorted spiral with the radius 
that decreases together with the temperature growth. 
One can also see, that the effective electron-electron in- 
teraction is attractive (Re[A (u)] > 0) in the range of the 
frequencies from zero to ^ O.SSrJmax- 

The free energy difference between the super- 
conducting and normal state (AF) can be calcu- 
lated on the basis of the formula [3]: 
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where p (0) denotes the value of the electronic density of 
states at the Fermi energy. The symbols Z,^ and rep- 
resent the wave function renormalization factors for the 
superconducting (S) and normal (N) state respectively. 

With the help of AF the specific heat in the super- 
conducting (C^) and normal (C^) state, as well as, the 
thermodynamic critical field (He) have been determined 
f3|. In Fig. [5] (A)-(B) we have presented the depen- 
dence of the specific heats and the thermodynamic crit- 
ical field on the temperature. On the basis of deter- 
mined thermodynamic functions the values of the ratios 
i?2 = AC (Tc) /C^ (Tc) and R3 = TcC" (Tc) /Hi (0) 
have been calculated. We have obtained: R2 = 1.65 and 
i?3 = 0.16. We notice, that the BCS model predicts: 
[^2]bcs ~ ^■^'^ ^"^^ [^sIbcs ~ 0.168 0. It is easy to see, 
that for Sulphur, R2 and R^ differ from the BCS values. 

Finally, we have calculated the ratio of the electron 
effective mass (m*) to the bare electron mass (mg): 
TO*/me = Re[Z(0)]. It has been stated, that m* is 
relatively high in the full range of the considered tem- 
peratures and ["^e/'^elmax = l-'''7 for T — Tc- 



IV. SUMMARY 

The high-pressure superconducting state in Sulphur 
have been analyzed in the framework of the Eliashberg 
approach. It has been stated that the exact values of 
the thermodynamic parameters cannot be calculated by 
using of the simple BCS model. 
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